Non-magnetic impurities to induce magnetism in a-PbO crystal structure 
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A new route to d° magnetism is established with help of the first principles methods. Non- 
magnetic interstitial impurities from group 14 in the periodic table are found to induce p-orbital 
magnetism in polycrystalline PbO-type structures. The half-filled p-orbitals occupied by two elec- 
trons is generated on the impurity site for which the ferromagnetic state of high stability is guaran- 
teed by the first Hund's rule. Since the impurity is embedded between layers of the host, its atomic 
radius is a key to tune not only its solubility but also the magnetic behavior: the on-site stability 
of the spin polarized state grows with reduction of the atomic radius while losing in the long-rang 
order interactions. However, for impurities of smaller radius the weaker inter-site magnetic cou- 
pling can be compensated by their concentration as the impurity solubility limit is shifted to higher 
magnitudes. 
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The research on magnetic semiconductors has inten- 
sified in recent years due to requirements imposed by 
the rapidly developing field of spintronics [1]. Originally, 
magnetic semiconductors were created by doping conven- 
tional semiconductors with magnetic ions whose d or / 
orbitals are partially filled 0, yj . The enormous atten- 
tion given to the so-called " diluted magnetic semiconduc- 
tors" has been rewarded with the discovery of a mech- 
anism of the 'intrinsic' magnetism in semiconductors - 
defect-induced magnetism. Initially, in the semiconduc- 
tors doped with magnetic ions, the intrinsic defects were 
considered only to mediate the magnetic coupling be- 
tween localized spins occupying the partially filled d or / 
orbitals of ions thus contributing to the collective mag- 
netism effect [J-|6|. However, a better understanding of 
the defect properties has revealed that the defects gener- 
ate sp localized spins able to generate the magnetic phe- 
nomenon [thTH- The discovery of defect- induced mag- 
netism dubbed as d a magnetism, i.e. magnetism which 
occurs not due to partially filled d orbitals of the mag- 
netic ions, brought new impetus into field of magnetic 
semiconductors and more importantly in spintronics. 

In order to pursue the spintronics application, the mag- 
netic semiconductors are required to establish the mag- 
netic percolation persisting at the room temperature. 
Such percolation is governed by interaction of the lo- 
calized spins provided they show also the high on-site 
stability of the spin-polarized state. In practice to in- 
duce d° magnetism, the weaker are interactions between 
the localized spins, the higher concentration of the de- 
fects carrying the localized spins has to be adapted @, HI ■ 
However, raising a defect concentration is not always a 
straightforward solution Q. For those defects known 
to induce magnetism: the vacancy and substi- 

tutional defects fl3l - fT5j , the low limit to the defect con- 
centration even at the most favorable growth conditions 
is applied because of their rather high formation energy 
[T3 |. Moreover, defects especially in high concentration 
are not always mechanically tolerated by the crystal lat- 
tice of the host, not to mention that the defect-induced 
lattice perturbation may lead to unwanted changes in 
the electronic properties |3| ■ Therefore, although an idea 



of intrinsic magnetism in semiconductors is attractive, it 
requires some efforts to bring it to the level of practical 
applications. 

In this respect, our recent finding of new route to 
d° magnetism offers an elegant solution and promises a 
breakthrough in development of magnetic semiconduc- 
tor [l6j]. We found that the Pb interstitial defects in 
a-PbO are capable to induce d° magnetism [16]. The 
origin of magnetism upon bonding of the impurity with 
the host, is unique; the Pb interstitial of Pb.;:6s 2 6p 2 
outer shell utilizes its only Pb^:6s 2 electrons to connect 
to the host while leaving on site two unpaired 6p 2 elec- 
trons. The Hund's rule dictates the spin alignment of 
the 6p 2 electrons which manifests in on-site magnetic 
moment of 2.0 /iB and the high stability of the spin 
polarized state defined by difference in total energy of 
the anti- and ferromagnetic states (AFM and FM) as 
E M = E AFM ~ £ FM =0.235 eV [jj]). Thus, the Pb atom 
being nonmagnetic in nature gains magnetization upon 
embedding as the interstitial defect into the ct-PbO crys- 
tal lattice. In analogy with magnetic ions, magnetism 
occurs due to partially filled orbitals, but with one sub- 
stantial difference, here it is due to the p orbital. 

In terms of large lattice distortion generated by the 
defects in crystalline materials (only crystalline systems 
have been considered in all previous studies devoted to 
d° magnetism [Tol - tl5j ). the layered type of the crystal 
structure of cv-PbO helps to avoid this issue. The in- 
terstitial defects in the a-PbO lattice are accommodated 
within the interlayer space or on the surface of crystal- 
lites thereby minimizing the negative impact of defects 
on the crystal lattice of the host. It results in lower- 
ing of the formation energy of defects and promises a 
better mechanical tolerance of the host lattice to such 
defects. Moreover, the mechanism of bonding of the in- 
terstitial to the host is unique itself because only the s 2 
valence electrons from both the host and the defect sides 
are involved in the bonding. In particular, the Pbi:6s 2 
outer shell electrons of the interstitial through excitation 
to the Pbi^s 1 ^ 1 state participate in the bond formation 
with the host Pb:6s 2 electrons [l6|. Since in the band di- 
agram of the host [TtJ the bands near the gap edges are 
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FIG. 1: (a) The a-PbO crystal structure which contains the 
impurity interstitial atom, (b) The list of impurities used 
to generate magnetism (the atomic number on the top and 
the standard atomic weight at the bottom), (c) The Pb 
self-interstitial: the spin density map is plotted with iso- 
values of ±0.003 e/A 3 in Xcrysden for the energy range 
(Ed + -EV)±0.15 eV. It demonstrates formation of the FM 
state at the impurity site for which defect tails can be traced 
up to the last oxygen atom shown. 



formed mainly by the host Pb:6p and the 0:2p electrons, 
but doping involves only the s electrons it is almost non- 
invasive to the host band behavior other than appearance 
of the defect-induced bands. Therefore, taking all these 
into account, a lattice of the a-PbO-type is expected to 
tolerate the high concentration of Pb interstitials to be 
limited only by thermodynamics of the defect formation. 

Synopsis on the origin of the p-orbital magnetism in 
the cv-PbO crystal lattice: the lone pair of s 2 outer shell 
electrons from the host and the s 2 p 2 outer shell from 
the impurity are required. Because the presence of s 2 
outer shell electrons is a typical feature of a-PbO lat- 
tice, 'other' compounds of this type can be used as the 
host lattice. In fact, a-PbO crystal structure is well 
recognized in superconductivity [18|: Fe-pnistides (ba- 
sics are FeSe, FeAs), cuprate (basics are CuO, CuS), 
and lanthanum compounds (basics are LaF, LaO). How- 
ever, since for the defect-induced magnetism the wide 
gap semiconductor is preferable, the a-PbO compound 
is seems to be the best candidate anyway. In regards to 
impurity choice, the similar statement is applied; instead 
of Pb atom 'other' chemical elements possessing the s 2 p 2 



outer shell configuration can be used as the magnetic 
impurity. The selection of impurity can be even further 
extended to the chemical elements of the s 2 p 3 outer shell 
which were found to build the on-site magnetic moment 
of 1.0 /zB QJ]. However, closer look to s 2 p 3 systems has 
revealed the weak spin-polarization of the p 3 localized 
state (for the Sb impurity, the spin polarization energy is 
only 0.046 eV) that implies the low temperature limit on 
an appearance of magnetism making it less suitable for 
device applications. Therefore, in this work we consider 
the chemical elements of the s 2 p 2 configuration, i.e. from 
group 14th in the periodic table (see the list of elements 
in Fig. HJb)), as potential candidates to induce d° mag- 
netism in the a-PbO crystal lattice. Because we expect a 
difference in the atomic radius of these elements influenc- 
ing their magnetic behavior, the possibility to tune mag- 
netism is investigated here. If our hypothesis is proven 
true, this approach can open a wide perspective to design 
the desired magnetic strength in the a-PbO magnetic 
semiconductor. 

In our study we applied the generalized gradient ap- 
proximation (GGA) with the PBE parametrization [l9| 
provided by WIEN2k package for the density functional 
calculations [20| (augmented plane wave + local orbitals 
approach). The Pb:5p, 5d, 6s, 6p and 0:2s, 2p electrons 
have been treated as the valence electrons. The su- 
percell approach (RKmax=7) with sufficiently large su- 
percell of 108-atom size (3x3x3 array of the primitive 
unit cells) has been used for single impurity calculation 
while the 190-atom size (4x4x3) supercell for the in- 
teracting defects. For integration of the Brillouin-zone, 
the Monkhorst-Pack scheme of the 5x5x4 (or 4x4x2) 
k-mesh was applied. The formation energies of the in- 
terstitial defects to be considered here for the vacuum 
conditions have been calculated with help of formalism 
described elsewhere [2l|. Because the distance between 
layers in the cv-PbO structure defined by the lattice con- 
stant Co contributes into the formation energy, the experi- 
mental value of co=5.07A [12] has been used for those cal- 
culations, while the band structure properties have been 
evaluated with the optimized lattice parameters as de- 
scribed in [2l[ . 

Our previous study [l6| has revealed that magnetism 
originated from the p-type localized spins of the Pb in- 
terstitial seems to combine advantages of the vacancy 
induced d° magnetism @-[ll| and magnetism due to the 
magnetic ions [1, Q . Because spin ordering of Qp 2 +y elec- 
trons on site of the Pb interstitial is governed by the 
Hund's rule, it therefore causes the high on-site stability 
of the ferromagnetic state. At the same time, hybridiza- 
tion of impurity state with the host lattice results in the 
extended defect tails promising to induce the long-range 
order interactions. As shown in Fig.[TJc), the defect tails 
appear in upper and lower layers, they are extended up 
to seven nearest-neighbors and show the higher spin lo- 
calization at the oxygen atoms. Although the impurity 
interacts with the top layer through bonding while with 
bottom layer only through the hybridization interactions, 



3 



TABLE I: The stability of the spin-polarized state as a func- 
tion of impurity atomic radius Ri m [2J| : the spin-polarization 
energy Em = -Eafm — Efm and the energy splitting of the 1" 
and l d bands (Ei — £2). Im-Pb is the double bond length. 



spin-up spin-down 



Im 


Rim, A 


Im-Pb, A 


E M , eV 


(E 1 -E 2 ), eV 


Pb 


1.81 


2.90 


0.235 


0.523 


Sn 


1.72 


2.86 


0.258 


0.585 


Ge 


1.52 


2.70 


0.306 


0.680 


Si 


1.46 


2.65 


0.338 


0.734 


c 


0.90 


2.30 


0.538 


1.058 



the defect tail are observed to be more pronounced at 
the bottom layer. It occurs because the Pb interstitial 
is tightly sandwiched between layers that results in its 
strong hybridization with the bottom layer. The redis- 
tribution of the spin density from defect site to the host 
lattice is probably reason of the on-site stability of the fer- 
romagnetic state found for the Pb interstitial (i?M=0.235 
eV) to be at least twice lower than that known for the 
magnetic impurity of d or / types [23] . It is natural to 
assume that for impurity with smaller atomic radius the 
hybridization with the bottom layer will be suppressed 
that should result in enhancement of on-site stability of 
the localized spins. Our study on magnetic behavior of 
the chemical elements from the group 14th embedded as 
the interstitial defects into the a-PbO crystal structure 
supports this hypothesis. We have considered an alter- 
ation in the band behavior as shown in Fig. [2] while the 
energetic (spin-polarization energy) and structural pa- 
rameters are disclosed in Table Q] 

In the band diagram (see Fig. [2] (a)), there are sev- 
eral bands associated with impurity, 1" and l d are those 
induced by p- localized electrons (the spin- up band 1" is 
occupied by two electrons, while spin-down band l d is 
empty) while the 2 U and 2 d bands forming the top of the 
conduction band (Ec) both appear due to the antibond- 
ing orbitals of the impurity- host bond (Im-Pb). It was 
found that reduction in atomic radius of dopant leads to 
shift of both impurity bands l u and l d (l u is occupied 
by two electrons, while l d is empty) towards the valence 
band (Ey) but in the same time the gradual enhance- 
ment of the energy splitting of these bands occurs. The 
2 U band is also shifted towards the lower energy due to 
shortening of the Im-Pb bond, while the energetic posi- 
tion of the 2 d band deviates in narrow energy range near 
~ 1.8 eV. The enhanced splitting of the l u and l d bands 
(see {Ei — E 2 ) in Table Ql indicates a gain in stability of 
the spin-polarized state and, indeed, a raise in the spin- 
polarization energy Em is observed. We found that for 
the Si and C impurities the on-site stability already ap- 
proaches the values known for the magnetic ions of d or 
/ types H|. 

In accordance with our expectation, a rise in on-site 
spin stability with reduction of the atomic radius of impu- 
rity is obtained. Moreover, the impurity of smaller radius 
was found to induce less distortion into the a-PbO crys- 
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FIG. 2: (a) The band diagram for the a-PbO crystal struc- 
ture containing the Pb self-interstitial defect: V and l d are 
the bands formed by the p-localized electrons Pb:6p^ +a while 
2" and 2 d are the antibonding orbitals of the Im-Pb bond 
0. (b) The energetic location of the 1", l d , 2" and 2 d 
bands relative the top of the valence band Ey for the differ- 
ent impurities. 



tal structure: in order to accommodate the Pb interstitial 
the layers of the host move apart while for impurity of 
smaller radius the interlayer distance is preserved. As a 
result, reduction in atomic radius of the interstitial atom 
results in suppression of the formation energy of defect: 
1.23 eV is found for Pb interstitial, 0.79 eV for the Sn 
interstitial, just above zero for Ge and Si, and it is becom- 
ing negative for C interstitial. The low formation energy 
promises not only the better mechanical tolerance of the 
lattice to defects but also the higher finite defect concen- 
tration to be reached: thermodynamically granted defect 
concentration for the negative formation energy can be as 
high as number of the sites available for bonding ~10 22 
cm~ 3 . Therefore, the C interstitial shows the most en- 
couraging magnetic behavior so far: the highest on-site 
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stability (Em =0.538 eV) and the negative formation en- 
ergy. Unfortunately, the high on-site stability suggests 
the stronger localization of the spin density at the defect 
site that might affect the long-order interactions. 

The collective magnetic ordering can only occur when 
two impurities are close enough to each other to establish 
the magnetic coupling of their localized spins and in this 
respect the long range order interactions play the cru- 
cial role. The magnetic coupling between impurities has 
been simulated for the system containing two interstitials 
of the same type placed on distance of ~4 A within the 
3x3x3 supercell. The p 2 +y electrons localized on impu- 
rities have been aligned on-site, while their inter-impurity 
ordering has been switched from antiferromagnetic to fer- 
romagnetic in order to evaluate Em — Eafm — £fm- The 
6Px +y state is exactly half filled and, therefore, for two 
interacting impurities when their localized electrons are 
ferromagnetically coupled (in this particular case all four 
localized spins in the system are aligned and the mag- 
netic moments is 4 /Ltf,), the inter-site virtual hopping is 
not allowed [!, [l(| ■ In case of AFM coupling, the virtual 
hopping is allowed leading to lowering of the total energy 
of the system such as the AFM state becomes the ground 
state. For two Pb interstitials placed on distance ~4 A we 
found that i?M(Pb-Pb)=-0.96 eV (a negative sign indicates 
the AFM ground state) while for the case of two interact- 
ing C interstitials it is only Bm(C-C) = "0-38 eV. The elec- 
tronic interactions and so the magnetic coupling between 
defects exponentially drop down as the distance between 
them increases. Thus, for two defects separated by a 
distance 11.5 A (for this calculation the size of the su- 
percell was 4x4x3), the energy difference between AFM 
and FM states is drastically supressed to i?M(Pb-Pb)=- 
0.0056 eV and E U (c-c) =-0.0023 eV for the Pb and C 
interstitials, respectively. Therefore, we found that the 
larger is the atomic radius of impurity, the stronger is the 
hybridization of the impurity state with the host lattice 
that results in extended defect tails being an essential 
to establish the long-range order interactions. For the 
magnetic percolation to be reached, a reduction in inter- 
impurity coupling to occur for the impurities of smaller 
atomic radius can be compensated by their concentration 
as the thermodynamic limit to the defect concentration 
is lifted to the higher magnitudes. 



In summary, we propose to generate d° magnetism in 
compounds of a-PbO crystal structure by doping with 
non-magnetic impurities belonging to group 14th of the 
periodic table. In analogy with magnetic ions, origin of 
magnetism is due to partially filled orbitals, but instead 
of d or / types it is being of p-type in this case. The 
partial occupation of p orbital appears as a result of the 
unique bonding of non-magnetic impurities of s 2 p 2 outer 
shell with the host lattice: two p 2 +y electrons are left 
on-site of impurity after bonding and the first Hund's 
rule dictates their ferromagnetic ordering. We found that 
the magnetic behavior of the dopants depends on their 
atomic radius: dopant with the smaller radius establishes 
the higher on-site stability of the localized spin. For the 
Si and C impurities, the on-site stability reaches such a 
high magnitude that it becomes c omp arable with that for 
the magnetic ions of d or / types [23j . Another benefit of 
the small atomic radius of impurity is its near zero defect 
formation energy that shifts the solubility limit to mag- 
nitudes as high as ~10 22 cm~ 3 . However, in terms of the 
collective magnetic ordering to occur, the trade-off be- 
tween the on-site stability and the long-range order inter- 
actions has to be made: the stronger on-site stability im- 
plies the weaker long-range order interactions. Therefore, 
to establish the magnetic percolation with non-magnetic 
impurities from 14th group of the periodic table, one rule 
should be followed: the smaller the atomic radius of el- 
ement the weaker magnetic coupling between impurities 
would occur to be compensated by higher doping concen- 
tration. 
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